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r~ i Abstract 

Nonlocal Lotka-Volterra models have the property that solutions concentrate as Dirac 
masses in the limit of small diffusion. Is it possible to describe the dynamics of the lim- 
iting concentration points and of the weights of the Dirac masses? What is the long time 
asymptotics of these Dirac masses? Can several Dirac masses co-exist? 

We will explain how these questions relate to the so-called " constrained Hamilton- Jacobi 
equation" and how a form of canonical equation can be established. This equation has been 
established assuming smoothness. Here we build a framework where smooth solutions exist 
J> ■ and thus the full theory can be developed rigorously. We also show that our form of canonical 

OO | equation comes with a structure of gradient flow. 

Numerical simulations show that the trajectories can exhibit unexpected dynamics well 
explained by this equation. 

Our motivation comes from population adaptive evolution a branch of mathematical 
ecology which models darwinian evolution. 



1 Motivation 

The nonlocal Lotka-Volterra parabolic equations arise in several areas such as ecology, adaptative 
dynamics and as limits of stochastic individual based models in the limit of infinite population. 
The simplest example assumes global interactions between individuals with a trait x and reads 

d t n e - eAn e = ^-R(x, I e (t)) , t > 0, x £ R d , (1.1) 
with a nonlinearity driven by the integral term 

J e (t) = / i/i(x)n e (t,x)dx. (1.2) 
Another and more interesting example is with local interactions 

d t n e (t,x) = -n e (t,x) ( r(x) - [ C (x , y)n e (t , y) dy] + eAn £ (t,x). (1.3) 
e V ./rod / 
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1 MOTIVATION 



We denote by n° > the initial data. 

These are called 'mutation-competition' models because the Laplace term is used for model- 
ing mutations in the population. Competition is taken into account in the second model by the 
competition kernel C(x,y) > and in the first model by saying that R can be negative for I t 
large enough (it a measure of how the total population influences birth and death rates). Such 
models can be derived from stochastic individual based models in the limit of large populations, 
[9j[6j[7]. There is a large literature on the subject, in terms of modeling and analysis, we just 
refer the interested reader to [12l 13 [23j 12S] . 

We have already normalized the model with a small positive parameter e since it is our goal 
to study the behaviour of the solution as e — > 0. The interesting qualitative outcome is that 
solutions concentrate as Dirac masses 

n € (t, x) « p(t)8(x — x(t)) . 

For equation (|l.ip . we can give an intuitive explanation; in this limit, we expect that the relation 
n(t, x)R(x, J(i)) = holds. In dimension 1 and for x i— > R(x,I) monotonic, there is a single 
point x = X(I) where R will not vanish and, consequently, where n will not vanish. A priori 
control of the total mass on n from below implies the result with x(t) = X(I(t)). 

In several studies, we have established these singular limits with weak assumptions 124] . 
A main new concept arises in this limit, the constrained Hamilton- J acobi equation introduced 
in |13| which occurs by some kind of real phase WKB ansatz (as for fronts propagations in 

nana nil) 

n e (t,x) = e u ^' x)/£ . (1.4) 
Here we have in mind the simple example of Dirac masses approximated by gaussians 

x-x\ 2 /2e _ e (-\x-x\ 2 -tln(2-we))/2t 

It is much easier to describe the limit of — | x — x\ 2 — eln(2-Ke)\ Dirac concentration points are un- 
derstood as maximum points of u e (t, x) in (jl.4j) . As it is well understood, these Hamilton-Jacobi 
equations develop singularities in finite time [TJ [JJl [TP] which is a major technical difficulty both 
for proving the limit and for analyzing properties of the concentration points x(t). 

This method in [13] of using eln(n e ) to prove concentration has been followed in several sub- 
sequent studies. For long time asymptotics (and not e — > but the two issues are connected as 
we explain in section [2|) it was introduced in [10] and used in \27\ 126] . More recently in [8] the 
authors come back on the Hamilton-Jacobi equation and prove that it makes sense still with 
weak assumptions for several nonlocal quantities Ik = f t^kn e (t, x)dx which can be characterized 
in the limit. 

Here we take the counterpart and develop a framework where we can prove smoothness of the 
various quantities arising in the theory. This opens up the possibility to address many questions 
that seem impossible to attain directly 

• Do the Dirac concentrations points appear spontaneously at their optimal location or do they 
move regularly? 

• In the later case, is there a differential equation on the concentration point x(t)7 It follows 
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2 A SIMPLE EXAMPLE: NO MUTATIONS 



from regularity that we can establish a form of the so-called canonical equation in the language 
of adaptive dynamics |12[ [TT] . This equation has been established assuming smoothness in [13] , 
in our framework it holds true. 

• In higher dimensions, why is a single Dirac mass naturally sustained (and not the hypersurface 
R(-,I) = for instance)? The canonical equation enforces constraints on the dynamics which 
give the explanation. 

• What is the long time behaviour of the concentration points x{t)l A simple route is that the 
canonical equation comes with the structure of a gradient flow. 

We develop the theory separately for the simpler case of global interactions model (jl.ip and 
or the global interaction model (|1.3p . For global interactions we rely on assumptions stated in 
section [3] and we give all the details of the proofs in the three subsequent sections. We illustrate 
the results with numerical simulations that are presented in section We give several extensions 
afterwards; in section [8] we treat the case with non-constant diffusion, and finally the case of 
local interactions in section [9j 

2 A simple example: no mutations 

The Laplace term in the asymptotic analysis of (jl.ip and (jl.3p is at the origin of several as- 
sumptions and technicalities. In order to explain our analysis in a simpler framework, we begin 
with the case of the two equations without mutations set for t > 0, x E 

d t n = nR(x,I(t)), I(t) = tp(x)n(t,x)dx. (2.1) 

d t n{t, x) = n(t, x) (r{x) - C{x, y)n(t, y) dy^j := n(t, x)R(x, I(t, x)) . (2.2) 

Also, we give a formal analysis, that shows the main ideas and avoids writing a list of assump- 
tions; those of the section [3] and [9] are enough for our purpose. 

In both cases one can easily see the situation of interest for us. The models admit a continuous 
family of singular, Dirac masses, steady states parametrized by y E M. d and the question is to 
study their stability and, when unstable, how the dynamics can generate a moving Dirac mass. 

The Dirac steady states are given by 

n(x;y) = p(y)5(x - y). 
The total population size p(y) is defined in both models by the constraint 

R{y,L(y))=0, 

with respectively for (|2.ip and (|2.2p 

L(y)) =i>{y)p(y), resp. l(y) = p{y)C(y,y). 

A monotonicity assumption in / for model (|2.ip . namely Ri(x,I) < shows uniqueness of I_(y) 
for a given y. In case of (|2.2p it is necessary that r(y) > for the positivity of p(y). 
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In both models a 'strong' perturbation in measures is stable, i.e. only on the weight; for 
n° = p°5(x — y), the solution is obviously n(t,x) = p(t)5(x — y) with 

-p(t) = p(t)R(y,il>(y)p(t)), resp. -p(t) = p(t)[r(y) - p(t)C(y,y)], 

and 

p(t) —> p(y). 

t— >oo — 

This simple remark explains why, giving I now, the hypersurface {x, R(x, I) = 0} is a natural 
candidate for the location of a possible Dirac curve (see the introduction). 

Apart from this stable one dimensional manifold, the Dirac steady states are usually unstable 
by perturbation in the weak topology. A way to quantify this instability is to follow the lines of 
|13j and consider at t = an exponentially concentrated initial data 

n°(x) = e u ° {x)/e — > p(x°)5(x - x°). 

It is convenient to restrict our attention to u® uniformly concave, having in mind the gaussian 
case mentioned in the introduction. Then, e measures the deviation from the initial Dirac state 
and to see motion it is necessary to consider long times as t/e or, equivalently, rescale the 
equation as 

e d t n e = n e R(x,I e (t,x)), 

and our goal is to prove that 

nJt, x) — > nit, x) = p(t)S(x — x(t)) . 

Also the deviation to a Dirac state turns out to stay at the same size for all times and we can 
better analyze this phenomena using the WKB ansatz (|1.4|) . Indeed, u e satisfies the equation 

d t u e = R(x,I e (t,xj). 

As used by [10], because u° is concave, assuming x h-» R(x,I e (t,x)) is also concave (this only 
relies on assumptions on the data), we conclude that u e (t, •) is also concave and thus has a unique 
maximum point x t (t). The Laplace formula shows that, with x(t) the strong limit of x e (t), 

n e (tx) _^ s , x _ 2 ( t) y 
J md n e (t,x)dx e^o 

With some functional analysis, we are able to pass to the strong limit in I £ and u t . Despite 
its nonlinearity, we find the same limiting equation, 

d t u = R(x, I(t, x)) , u(t = 0) = u°. (2.3) 

Still following the idea introduced in [13], we may see I(t, x) or pit) as a Lagrange multiplier for 
the constraint 

maxu(t,x) = = u(t,x(t)), (2.4) 
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3 GLOBAL COMPETITION: ASSUMPTIONS AND MAIN RESULTS 



which follows from the a priori bound < pit) < pu < oo. The mathematical justification of 
these developments is rather easy here. For model (I2.ip it uses a BV estimate proved in [4]. 
For model (I2.2D one has to justify persistence (that is p e stays uniformly positive) and strong 
convergence of p e {t). All this work is detailed below with the additional Laplace terms. 

The constraint (12. 4p allows us to recover the 'fast' dynamics of I{t) and p(t). Indeed, combined 
with ([OD , it yields 

R(x(t),I(t)) =0, resp. R(x(t), I(t,x(t))) = 0. (2.5) 

Assuming regularity on the data, u(t, x) is three times differentiable and the constraint (12, 4p 
also gives 

Vu(t,x(t)) = 0. 

Differentiating in t, we establish the analogue of the canonical equations in [TTJ (see also [5j 1131 

E3ETJ) 

2{t) = (-D 2 u(t,x(t))y 1 .D x R(x(t),I(t,x(t))), x(0) = x°, (2.6) 

(in the case of model (12, 2p . this means the derivative with respect to x in both places). Inverting 
/ from the identities (|2,5p gives an autonomous equation. 

This differential equation has the structure of a gradient flow and this makes easy to analyze 
its long time behaviour, that is equivalent to know what are the stable Dirac states for the weak 
topology; the so-called Evolutionary attractor or Convergence Stable Strategy in the language 
of adaptive dynamics [12\ 127] . Eventhough this is less visible, it also carries regularity on the 
Lagrange multiplier p(t) which helps for the functional analytic work in the case with muta- 
tions. Another use of (|2.6p is to explain why, generically, only one pointwise Dirac mass can be 
sustained; as we explained earlier, the equation on x(t) also gives the global unknown I(t) by 
coupling with (|2.5p and this constraint is very strong. See section [7] for an example. 

To conclude this quick presentation, we notice that the time scale (in e here) has to be precisely 
adapted to the specific initial state under consideration. The initial state itself also has to be 
'exponentially' concentrated along with our construction; this is the only way to observe the 
regular motion of the Dirac concentration point. This is certainly implicitly used in several 
works where such a behaviour is displayed, at least numerically. Of course, there are many 
other ways to concentrate the initial state with a 'tail' covering the full space so as to allow that 
any trait x can emerge; these are not covered by the present analysis. 

3 Global competition: assumptions and main results 

As used by |10j , concavity assumptions on the function u t in (jl.4p are enough to ensure concen- 
tration of n e (t, ■) as a single Dirac mass. We follow this line and make the necessary assumptions. 
We start with assumptions on ip: 

o < < V> < V>m <oo, i; ew 2 >°°(m. d ). (3.1) 

The assumptions on R £ C 2 are that there is a constant Im > such that (fixing the origin 
in x appropriately) 

max R(x,Im) = = R(0,I M ), (3.2) 

x€M. d 
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-K^ 2 <R{x,I)<K -K 1 \x\ 2 , for < I < I M , (3.3) 

- 2K X < D 2 R(x,I) < -2Ki < as symmetric matrices for < I < I M , (3.4) 

-K 2 < — < -K 2 , A(^R) > -K 3 . (3.5) 
At some point we will also need that (uniformly in < I < Im) 

D 3 R(-,I) € L°°(R d ). (3.6) 

Next the initial data n® has to be chosen compatible with the assumptions on R and if). We 
require that there is a constant 1° such that 

< 1° < / e (0) := f ifi(x)n°(x)dx < I M , (3.7) 

that we can write 

n ° = e u ° /e , with u° e G C 2 (M d ) (uniformly in e), 

and we assume uniform concavity on u e too. Namely, there are positive constants L , Lq,L 1: L\ 
such that 

-La-Li\A 2 < u° € (x) < L -Li\x\ 2 , (3.8) 

- 2L t < D 2 u° e < -2L~i. (3.9) 



For Theorems 13.21 and 13.31 we also need that 

D 3 u® G L°°(M. d ) componentwise uniformly in e, (3.10) 
n®(x) — > p 5{x — x ) weakly in the sense of measures. (3-11) 

Next we need to restrict the class of initial data to fit with R through some compatibility 
conditions 

AL 2 X < Ki < K x < 4L_l. (3.12) 



In the concavity framework of these assumptions, we are going to prove the following 

Theorem 3.1 (Convergence) Assume (^Tj) - (|53|) . (f377|) - (f3T9|) and (f37T2|) . Then for allT > 0, 
there is a eo > such that for e < eo an d t G [0, T], i/ie solution n e to (jl.ip satisfies, 

0<p m < Peif) ■= / n e < p M + Ce 2 , < I m < I e (t) < I M + Ce 2 a.e. (3.13) 

for some constant p m , I m . Moreover, I e is uniformly bounded in BV(M + ) and after extraction 
of a subsequence I t 

J e (t)— )-J(t) in L} oc (R + ), I m <I(t)<I M a.e., (3.14) 

e— >0 
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and I(t) is non- decreasing. We also have weakly in the sense of measures for a subsequence n t 

n e (t,x) — > p(t) 6(x-x(t)). (3.15) 

Finally, the pair (x(£),/(t)) also satisfies 

R(x(t),I(t)) = a.e. (3.16) 

In particular, there can be an initial layer on I e that makes a possible rapid variation of I e at 
t ~ so that the limit satisfies R(x° , I(0 + )) = 0, a relation that might not hold true, even with 
O(e), at the level of n t . 

Theorem 3.2 (Form of canonical equation) Assume (|3.ip - (|3.12p . Then, x(-) belongs to 
W 1>oc (R + ;R d ) and satisfies 

x(t) = (-D 2 u(t,x(t)))~ 1 .V x R(x(t),I(t)), x(0) = x°, (3.17) 

with u(t,x) a C 2 '-function given below in (|5.8p . D 3 u G L oc (W 1 ), and initial data x° given in 
(l3~TTjl . Furthermore, we have I(t) G VF 1 '°°(IR + ). 

We insist that the Lipschitz continuity at t = is with the value 1(0) = lim t ^ + I(t) ^ 
lim^o-^; the equality might hold if the initial data is well-prepared. 

Theorem 3.3 (Long-time behaviour) With the assumptions ()3. lj) - (j3. 12|) . equation (|3. 17[) 
has a structure of gradient flow, the limit I(t) is increasing and 

I(t) — > I M , x(t) — »• xoo = 0. (3.18) 

t— >oo t— >oo 
Finally, the limit is identified by Vi?(xoo = 0,1m) = (according to (13.2 ). 

It is an open question to know if the full sequence converges. This is to say if the solution 
to the Hamilton-Jacobi equation is unique. The only uniqueness case in [1] assumes a very 
particular form of R(-, ■). 



4 A-priori bounds on p £ , I £ and their limits 

Here, we establish the first statements of Theorem 13.11 As in [3] we can show with (|3.2p and 
(|331) that I e < I M + Ce 2 . With (j3TTT) (the bounds on V), we also have that 

Pe(t) <hi/ipm + Ce 2 . (4.1) 

To achieve the lower bound away from is more difficult. We multiply the equation (jl.ip by 
ip and integrate over R rf , to arrive at 



^j-Ie{t) = \ J ^ Rn £ dx + e J n e Aip dx. (4.2) 
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1 f 

We define J e (t) := — t/jRn e dx and calculate its time derivative 

— J e (t) = - J ipR (—Rn e + eAn^j dx H — J ^;n e -^y (^J e + e J n e Aipdy^J dx. 
So we estimate it from below using (13. If) , (13.5f) and (I4.ip by 

^ J e (f) > -C + i J e (t) J i>nj^ dx, 
and we may bound the negative part of J e by 

^{Ut))-<C-—h{t){J e {t))-- (4-3) 
dt e 

Now for e small enough, we can estimate I e (t) as 

I e (t) = J £ (0) + T 7 e (s) = 7 e (0) + f J t {s) ds + 0(e) > 7°/2 - f\j e (s))- da, (4.4) 
and plugging this in the estimate (I4.3P leads to 



J t {Je{t))- <C-^ ll° /2 - j| (J e (,))_ eta 1 (J e (t))_. 

Now for T > fixed. If there exists T' <T such that j Q (J e (s))_ (is = 7 /4, then we have 

^(J e (t))_ <C- ^(J e (t))_, <t <T'. 

Thus we obtain 

(J e (t))- < (Mt = 0))_e-^°W + J^L (l - e-^/W) . 

Then for e < £o{T) small enough, we conclude that such a T' does not exist i.e. 

[ (J £ ( S ))_d S <7°/4. (4.5) 
Jo 

So from (14. 4|) . we obtain 

7 e (i) > 7°/4, (J e (*))- ^0 a. e. in [0,T]. (4.6) 

e->0 

This also gives the lower bound p m < p e (t) with p m := 7°/(4V>m)- 

Finally, the estimate (|4.5p and the 7°° bounds on 7 e (t) give us a local BV bound, which will 
eventually allow us to extract a convergent subsequence for which (13. 14f) holds. The obtained 
limit function 7(t) is non-decreasing because in the limit the right-hand side of (|4.2p is almost 
everywhere non-negative. 
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5 Estimates on u e and its limit u 

In this section we introduce the major ingredient in our study, the function u € := eln(n e ). We 
calculate 

dtn e = n e dtu e /e, Vn e = n e Vu e /e, An e = n e Au e /e + n e \Vu e \ 2 /e 2 . 
Plugging this in (jl.ip . we obtain that u e satisfies the Hamilton- Jacobi equation 
d t u e = \Vu e \ 2 +R(x,I e (t)) + eAu e , x E M d , t>0, 

(5.1) 

u e (t = 0) = eln(n°) := u°. 

Our study of the concentration effect relies mainly on the asymptotic analysis of the family 
u e and in particular on its uniform regularity. We will pass to the (classical) limit in (|5.ip . and 
this relies on the 

Lemma 5.1 With the assumptions of Theorem \3.1[ we have for t > 0, 

-Lo ~Li\x\ 2 - &dL x t < u € (t,x) < L - Ti\x\ 2 + (K + 2deLi) t, (5.2) 

-1L X < D 2 u e (t,x) < -2L~i. (5.3) 

This Lemma relies on a welknown (and widely used) fact that the Hamilton-Jacobi equations 
have a regime of regular solutions with concavity assumptions, [U [19] . 

5.1 Quadratic estimates on u e 

First we achieve an upper bound, defining u e (t,x) := Lq — L\\x\ 2 + Co(e)t with Cb(e) := Kq + 
2deL u we obtain thanks to (pT3l) . (J3T81) and (I3TT21 that u e (t = 0) > u° e and 

9 t H e - |Vn e | 2 - R(x,I e ) - eAu e > C (e) - 4Ll\x\ 2 -K + K 1 \x\ 2 - 2deL 1 > 0. 

Next for the lower bound, we define u e (t, x) := — L — L_i\x\ 2 — eC\t with C\ := 2dL 1 , we have 
u^(t = 0) < u° e and 

d t u e - \Vu e \ 2 - R(x,I e ) - eAu e < -ed - 4L?|x| 2 + K^xf + e2dL l < 0. 

This concludes the proof of the first part of Lemma 15. II i.e. inequality 15.21 

5.2 Bounds on D 2 u e 

We show that the semi-convexity and the concavity of the initial data is preserved by equation 
(|5.ip . For a unit vector £, we use the notation := V^u e and := V 2 ^u e to obtain 

u^t = R^(x, I) + 2Vu ■ Vu^ + eAu^, 

v,££t = R^(x, I) + 2Vii£ ■ Vii£ + 2Vu • Vu^g + eAu££. 
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The first step is to obtain a lower bound on the second derivative i.e. semi-convexity. It can 
be obtained in the same way as in [24J: Using \S/u^\ > \u^\ and the definition w(t,x) := 
min^u^(t, x) leads to the inequality 

dm > -2K X + 2w 2 + 2Vu ■ Vw + eAw. 

By a comparison principle and assumptions (|3.9|) . (|3,12p . we obtain 

w>-2L v (5.4) 

At every point (t, x) £ R + x M. d , we can choose an orthonormal basis such that D 2 u e (t,x) is 
diagonal because every symmetric matrix can be diagonalized by an orthogonal matrix. So we 
can estimate the mixed second derivatives in terms of u^. Especially we have |Vn^| = \u^\. 

This enables us to show concavity in the next step. We start from the definition w(t, x) := 
max£it££(t,x) and the inequality 

d t w < -2K 1 + 2w 2 + 2Vu ■ Vw + eAw. 

By a comparison principle and assumptions f)3.9f) . (13. 12H . we obtain 

w<-2L l . (5.5) 

From the space regularity gained and (|3.9p . we obtain Vu e locally uniformly bounded and 
thus from (15. 1|) for e < eo that &tu e is locally uniformly bounded. 



5.3 Passing to the limit 

From the regularity obtained in section 15.21 it follows that we can extract a subsequence such 
that, for all T > 0, 

Ue (t,x) — > u(t,x) strongly in L°° (o,T;W^ (R d )^j , 



u e (t,x) — ^ u(t,x) weakly-* in L°° [0,T;W^(R d )) n W 1 ' 00 (0, T; L% c (R d ] 



and 

-La ~ Li\A 2 < u(t,x) < L - Li|x| 2 + K t, -2L X < D 2 u(t, x) < -2L X a.e. (5.6) 

ueW^(R + xR d ). (5.7) 

Notice that the uniform W 2 ^(R. d ) regularity also allows to differentiate the equation in time, 
and find 

= 0_ R fa 1{t )) + 2Vu.[VR(x, I(t)) + D 2 u.Vu]. 

This is not enough to have C 1 regularity on u. 
We also obtain that u satisfies in the viscosity sense (modified as in [HIM]) the equation 

§ i u = R[x,I{t)) + \Vu\ 2 , 

(5.8) 

max R d u(t, x) = 0. 
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The constraint that the maximum vanishes is achieved, as in [3], from the a priori bounds on I e 
and flggp . 

In particular u is strictly concave, therefore it has exactly one maximum. This proves (|3.15|) 
i.e. n stays monomorphic and characterizes the Dirac location by 

maxu(t,x) = = u(t,x(t)). (5.9) 

Moreover, as in [24] we can achieve (I3.16P at each Lebesgue point of I(t). 
This completes the proof of Theorem 13.11 

6 Canonical equation, time asymptotic 

With the additional assumptions ()3.6[) and (|3.10p . we can write our form of the canonical equa- 
tion and show long-time behavior. To do so, we first show that the third derivative is bounded. 
This allows us to establish rigorously the canonical equation while it was only formally given 
in |13|. 124] . Prom this equation, we obtain regularity on x and /. For long-time behavior we 
show that I is strictly increasing as long as Vi?(x, Tj ^ and this is based on the gradient flow 
structure of the equation. 

6.1 Bounds on third derivatives of u e 

For the unit vectors £ and rj, we use the notation := Vgu e , u^ v := V^ite and u^ v := V^ v u e 
to derive 

dtu^ = 4S7u^ n ■ Vug + 2Vu v ■ Vu££ + 2Vu • Vu^ v + + eAu^. 

Again we can can fix a point (t,x) and choose an orthonormal basis such that D 2 (V r? it e (t, x)) 
is diagonal. Let us define 

M\(t) := maxti^ r) (t, x). 
x,£ ,77 

Since — = V- V u^(t,x), we have M\(t) = max \u^ v (t, x)\. So with the maximum prin- 
ciple we obtain 

—Mi < 4iiMi||L> 2 u e || 00 + 2dM 1 \\D 2 u t \\ oc + R^. 

Assumption (I3.10p gives us a bound on M\{t = 0). So we obtain a L°°-bound on the third 
derivative uniform in e. 

6.2 Maximum points of u t 

Now we wish to establish the canonical equation. We denote the maximum point of u e (t, •) by 
x e (t). 

Since u e 6 C 2 , at maximum points we have Vu e (£, x e (t)) = and thus 

-Vu £ (t,x e (t)) =0. 
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6.3 Long-time behaviour 



6 CANONICAL EQUATION, TIME ASYMPTOTIC 



The chain rule gives 

Of 



Vu e (t,x e (t)) +Dlu e (t,x e {t))x e (t) = 0, 



and using the equation (|5.ip . it follows further that, for almost every t, 

d 

Dlu e (t,x e (t))x e (t) = -—Vu e (t,x e {t)) = -V x R(x e (t),I e {t)) -eAV, 



u f . 



Due to the uniform in e bound on D 3 u e and ReC 2 , we can pass to the limit in this equation 
and arrive at 



2(t) = (-D 2 u(x(t),t)) .V x R(x(t),I(t)) a. 



e. 



But we can obtain further regularity in the limit and obtain the equations in the classical 
sense. We first notice that, from R(x(t), T(t)) = and the assumptions (|3.3p . x(t) is bounded 
in So we obtain that x(t) is bounded in W 1,co (M. Jr ). Because / i— > R(-,I) is invertible, 

it follows that J(i) is bounded in W 1 '°°(R + ); more precisely we may differentiate the relation 
(I3.16P (because R G C 2 ) and find a differential equation on I(t) that will be used later: 

x-(t) ■ V X R + I(t) V/.R = 0. 
This completes the proof of Theorem I3.2L 

6.3 Long-time behaviour 

It remains to prove the long time behaviour stated in Theorem 13,31 
We start from the canonical equation 

j t x(t) = (-D 2 u)- 1 VR(x(t),I(t)), 
and use some kind of gradient flow structure. We calculate 

^R{x(t)J(t)) = VR(x(t), 1(f)) jx[t) + Rj (x(t),I(t)) ^ 



~dt' 



= VR(x(t), I(t)) (-D 2 U y 1 VR(x(t), I(t)) + Ri(x(t), I(t)) 
Now we also know from (|3.16p that the left hand side vanishes. Then, we obtain 

The inequality is strict as long as I(t) < Im- Consequently, we recover that I(t) is non- 
decreasing, as t — > oo, I(t) converges, and subsequences of x{t) converge also (recall that x{t) is 
bounded). But we discover that the only possible limits are such that Vi?(xoo, 1^) = 0. With 
relation (|3.16p . assumptions (|3.2[) and (|3.5j) this identifies the limit as announced in Theorem 
1331 



12 



7 NUMERICS 




(a) asymmetric I.D. (b) symmetric I.D. 

Figure 1: Dynamics of the density n with asymmetric initial data (|7.ip (left) and symmetric 
initial data (|7.2p (right). These computations illustrate the effect of the matrix (— D 2 u) 1 in 
the dynamics of the concentration point according to the form of canonical equation (|3.17p . 

7 Numerics 

The canonical equation is not self-contained because the effect of mutations appears through 
the matrix (— D 2 u)~ 1 . Nevertheless it can be used to explain several effects. The purpose of 
this section is firstly to illustrate how it acts on the dynamics, secondly to see the effect of e 
being not exactly zero, and thirdly to explain why it is generic, in high dimensions as well as in 
one dimension |24j that pointwise Dirac masses (and not on curves) can exist. 

We first illustrate the fact that a isotropic approximation of a Dirac mass will give rise to 
different dynamics than an anisotropic. This anisotropy is measured with u and we choose two 
initial data. In the first case —D 2 u° is "far" from the identity matrix and in the second case it 
is isotropic: 

no(x,y) = C mass exp(-(j; - .7) 2 /e - 5(y - .7) 2 /e), (7.1) 

n (x, y) = C m ass exp(-2.4(x - .7) 2 /e - 2A(y - .7) 2 /e). (7.2) 

We also choose a growth rate R with gradient along the diagonal: 

R(x,y,I)=2-I-(x + y). (7.3) 

Here although, we start with initial data centered on the diagonal and Vi? pointing along the 
diagonal to the origin, the concentration point leaves the diagonal with the anisotropic initial 
data (|7.ip (cf. Figure [Q (a)). The isotropic initial data moves along the diagonal as expected by 
symmetry reasons (cf. Figure Q] (b)). 

The numerics has been performed in Matlab with parameters as follows. The three plots cor- 
respond to t = 0, 40 and 80 in units of dt: e is chosen to be 0.005 and C mass such that the initial 
mass in the computational domain is equal to 0.3. The equation is solved by an implicit-explicit 
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finite-difference method on square grid consisting of 100 x 100 points (time step dt = 0.01). 

The second example is to illustrate the role of the parameter e for symmetric initial data: 

n {x,y) = C mass exp(-(x - .3) 2 /e - (y - .3) 2 /e), (7.4) 



R(x, y, I(t)) =0.9-1 + 5(y - .3) 2 + + 2.3(x - .3) with I(t) := J n(t, x) dx. 



(7.5) 



In this example, we start with symmetric initial data centered on the line y = 0.3 and the 
gradient of R along this line (y = 0.3) is (1,0). Hence, the canonical equation in the limit e = 
predicts a motion in the x direction on this line. One however observes in Figure [2] that the 
maximum point leaves this line because e does not vanish. Notice that R = below the line 
y = 0.3. 

In this computation, performed with Matlab, e is chosen to be 0.004, C mass such that the 
initial mass in the computational domain is equal to 0.3 and square grid consisting of 150 x 150 
points (time step dt = 8.8889 • 10" 4 ). 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
1 







30 



0.2 



Figure 2: This figure illustrates the effect of e being not exactly zero. The dynamics of the 
density n with symmetric initial data is plotted for t; = 0, 160 and 220 in units of dt and the 
limiting behavior is a motion along the axis y = 0.3 . 



With our third example we wish to illustrate that, except in particular symmetric geometries, 
only a single Dirac mass can be sustained by our Lotka-Volterra equations. We place initially 
two symmetric deltas on the x and the y-axis: 



no(x,y) = C n 



< xp ( -— ( (x - .25\/2) 2 + y 2 X) +exp (-— ( {y - .25\/2) 2 + 



x 



We seek for asymmetry in the growth rate R under the form 

R(x, y, I) = 3 - 1.5/ + 5.6(y 2 + R e x 2 ). 



(7.6) 
(7.7) 
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8 EXTENSION: NON-CONSTANT DIFFUSION 
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(a) asymmetric R (b) symmetric R 

Figure 3: This figure illustrates that, except particular symmetry conditions, a single Dirac 
mass is exhibited by Lotka-Volterra equations. We depict the density n with asymmetric (left) 
and symmetric (right) growth rate R plotted for t = 5, 90 and 180 in units of dt. 



In the special case R e = 1, all isolines of R are circles then the two concentration points just 
move symmetrically to the origin cf. [3] (b). However, if we choose R e = 1.1 i.e. all isolines 
of R are ellipses then one of the two concentration points disappears cf. Figure [3] (a). The 
intuition behind is that the canonical equation (I3.17h should hold for the two points. However 
the constraint (|3.16p given by p(t) is the same for the two points and this is a contradiction. 
One of the two points has to disappear right away. 

The numerics is performed with e = 0.003 and C mass such that the initial mass in the compu- 
tational domain is equal to 0.3. The equation is solved by an implicit-explicit finite-difference 
method on square grid consisting of 100 x 100 points (time step dt = 0.001). 

8 Extension: non-constant diffusion 

Our results can be extended to include a diffusion coefficient depending on x. This leads to the 
equation 

d t n t - eV • (&(jc) Vn e ) = ^R(x, I 6 (t)) , t > 0, x G R d . (8.1) 
Our assumptions on b are that there are positive constants b m , &Mj B%, B 2 and B% such that 

b m <b<b M , \Vb(x)\<B 1 -^—, \Tr(D 2 b(x))\ < B 2 ) \D 3 b\ < B 3 . 

1 + |x| 1 1 (1 + \x\y 1 1 

(8.2) 

Our assumptions on the initial data and on R are the same, as before (|3.ip - (|3.1ip . However, 
we have to supplement the assumption (|3.5p to take b into account: 

V6 • V(ipR) > -K 3 . (8.3) 

These assumptions will in the following allow us to obtain a gradient bound 

\Vu e (t,x)\<C Wu {l + \x\). (8.4) 
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This bound enables us to formulate the compatibility conditions which replace ()3.8[) : we need 

B 2 Cl u - 2K 1 < (8.5) 

and define 



_ 2B X - JAB\ - 2b M [B 2 C 2 u - 2Ki) 

K b : = 

OM 



to require 



-2B X - J ABf + 2b m {B 2 C 2 U + 2K X ) 
K b := 



-Kb<D z u" £ <-K b , (8.6) 
4b M L\ <Kx<K x < 4b m Lj. (8.7) 



Our goal is to prove the following 

Theorem 8.1 (Convergence) Assume flUD -QSM, (USD, (pT2|) . (lOjl . (|Q> . (18^1) and (IHTfj) . 

T/ien i/ie solution n e to (|8. 1 j) satisfies for all T > 0, /or e small enough and t G [0, T] 

< p m < p £ (t) < pjjf + Ce 2 , J m < I e (i) < I M + Ce 2 a.e. (8.8) 
Moreover, there is a subsequence I e such that 

I € (t)^I(t) inLj oc (R + ), I m <I(t)<I M a.e., (8.9) 

e— >U 

and J(i) is non- decreasing. Furthermore, we have weakly in the sense of measures for a subse- 
quence n e 

nJt,x) — > pit) 5(x - x(t)), (8.10) 

and the pair {x{t),I{t)j also satisfies 

R(x(t),I{t)) =0 a.e. (8.11) 

Theorem 8.2 (Form of canonical equation) With the assumptions (|3.ip - (|3.8p . (|3.10p . (|3.1ip . 
([83]) . (g^D , (USD, (US]) and ([821), xwo VK 1 '°°(R + )-/unc^on satisfying 

x(t) = (-D 2 u(t,x(t))y 1 .V x R(x(t),I(t)), x(0) = x°, (8.12) 

with u(t,x) a C 2 -function given below in ([QT|h D 3 u £ L°°(R d ), and initial data x° given in 
([3~TT]) . Furthermore, we have I(t) € W 1 '°°(R + ). 



The end of this section is devoted to the proof of these Theorems. The a priori bounds (j8.8|) . 
(|8.9p on p e and can be established as before. 

As before we study the function u e := e ln(n e ). We obtain that u e satisfies the Hamilton-Jacobi 
equation 

d t u £ = R(x,I e (t)) +b\X7u e \ 2 + e\7b- Vu e + ebAu e , t>0,x£R d , 
« e (t = 0) = eln(n°). [ j 

In order to adapt our method to this equation we need a bound on the gradient of u e . We 
achieve this following arguments in [20^ [3] : 
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gradient bound Let us define v(t,x) by u e = K v — v 2 where we choose K v large enough to 
have v > 5 > on [0,T] uniform in e. We calculate 

Vu e = -2vVv and Au e = -2vAv - 2\Vv\ 2 

and obtain from (|8. 13[) 

-2vd t v = R + 4£w 2 |Vt;| 2 - 2eV6 • vVv - 2evbAv - 2e6|Vi;| 2 . (8.14) 
Dividing by — 2v, taking the derivative with respect to X{ and defining p := Vi>, we have 

dm = ~ (j^j ~ 2pib\p\ 2 - 2vb Xi \p\ 2 - 4vbp ■ V Pi + eV^. • p + eVb ■ V Pi 

b Xi ,2 b i i2 b 
+ eb x Av + ebApi + e—\p\ - e—»\p\ Pi + 2e-p ■ Vpj. 

v v v 

Multiplying (18. 14H by — ^- and adding to the equation above, we obtain 
2bv 

b Xi v\ _ ( R\ „ ., |2 „ , , |2 



dt ( Pi - ) = - ( ^ ) - 2pi6|p| 2 " 2^M 2 " 4^6p • Vpi + eVft^ • p + eVb ■ Vp; 



+ e&Ay>j + e^|p| 2 - e\\p\ 2 Pi + 2e-p • Vj* + %^ + 26 !Ei v|p| 2 - e^Vb ■ p - e-^\p\ 2 . 
v v z v 2bv b v 

Now we define 

M b (t) := max [{ Pi )_, { Pi ) + } > 0. (8.15) 

l,X 

For maxj i2 ;(pj)_ < maxj iX (j9j)_|_, we have 

' ^ ^ ™ "^ 3 + 2|v||6 Xi |d 2 M b 2 

&aij .2,1^2 , n , oil ||„.| j2 ti^2 , JIT71.I )l>f. i J&Sil j2 7i>r2 



$ ( M 6 - ] < C - 2b m M b i + 2\v\\b Xt \d 2 M b z + ed|V6 X4 | M 6 

+ J-^d z M' b z + C + 2\b Xi \\v\d z M£ + e ^d\Vb\M b + e '-^d z M } 



b . v 

Since — — is bounded, we have M b bounded from above. 
b 

For maxj j2 .(pj)_ > max.i x (pi)+, we show similarly a bound on M b and therefore achieve ([87 
To prove the concavity and semi-convexity results, we only give formal arguments for the limit 
case. To adapt the argument for the e-case is purely technical: 
For a unit vector £, we define u% := V^u e and := V^u e to obtain 

d tUi = R^ + b^\Vu\ 2 + 2bVu-Vu{., (8.16) 

and 

dtUtf = R^ + b^\Vu\ 2 + 46 5 Vu • Vu ? + 2bVu ■ Vu^ + 26|Vn f | 2 . (8.17) 
With the definition w(t,x) := max^ u^(t,x) and assumptions (18.2j) we have 
a r W < -2Ki + J3 2 C^ U + 4Bi|uJ| + 2bVu ■ VuJ + 2b M w 2 . 
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With assumption (|8.5p . is a supersolution to 

d t w* = -2K 1 + B 2 C\ U - ±B x w* + 2bVu • Vw* + 2b M {w*) 2 , 
so we know from assumption (j8.6j) that uJ < 0. Therefore it follows further that 

w < Kb. 

For the lower bound, we use the definition w(t,x) := min^u^(t, x) and the inequality 
d t w > -2K_ X - B 2 C% U - AB x \w\ + 2bVu ■ Vw + 2b m w 2 . 
Since we already know that w < 0, we obtain 

w > K b . 

We can achieve this at the e-level using the equation 

d t Uj: = R% + 6 5 |Vn| 2 + 26V u ■ Vu ? + eVfe 5 • Vu + eVb ■ Vu ? + efe^Au + ebAu^, (8.18) 

and 

9ttt^ = % + &&|Vu| 2 + 46 5 Vn • Vu 5 + 26Vn • Vu^ + 26|Vu 5 | 2 

+ eVbtf ■ Vu + 2eV6 ? • Vu ? + eV6 • Vu^ + e& ?5 Au + 2e6 f Au f + ebAu^. (8.19) 



Now we define 



26 ? ?%-2& 2 
/ := — — and g :— 



multiply (I8.18j) by /, substract it from (18.19j) . multiply (18.13R by g, substract it to obtain 

d t(utf - fu£ - gu) = Rf-j. + 26V u ■ Vu^ + 26|Vu 5 | 2 + eVb^ ■ Vu 

+ 2eV6 ? Vu 5 + eVb ■ Vu^ + ebAu^ - fR$ - efVb^ ■ Vu - efVb ■ Vu ? - gR- egVb ■ Vu. 

(8.20) 

The remaining steps can be done similar as before. For the Hamilton-Jacobi-equation on u, we 
obtain the variant 

d t u = R(x,I(t)) + b(x) |Vu| 2 , 

maxu(i,x) = 0, Vt > 0. ^ 8 ' 21 ^ 



9 Local competition 

The other class of models we handle are populations with localized competition kernel C(x, y) > 
0, that is 

d t n e (t,x) = -n e (t,x) (r{x) - C(x,y)n e (t,y) dy) + eAn e (t,x). (9.1) 
e V JR d J 
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The term r(x) is the intra-specinc growth rate (and has a priori no sign) and the integral term 
models an additional contribution to the death rate due to competition between different traits. 
Notice that the choice C(x, y) = tp(y)^(x) will reduce this model to a particular case of those in 
(jl.ip . This class of model is also very standard, see [9j [22J [26] and the references therein. 

For the initial data, we assume as before (|3.8[) — (|3.1ip . Concerning r(x) and C(x, y) we assume 
C 1 regularity and that there are constants pu > 0, > 0... such that 

C(x,x)>0 Vx£R d , (9.2) 



/ / n(x)C(x,y)n(y)dydx > — [ n(x)dxj r(x)n(x) dx Vn G L+(R d ). (9.3) 

This assumption is weaker than easier conditions of the type 

C(x,y) > —r(x) or C(x,y) > - [r(x) +r(y)}. 

Pm *Pm 

Because it is restricted to positive functions, it is a pointwise positivity condition on C(x,y) in 
opposition to the positivity as operator that occurs for the entropy method in [18J . 

Then, we make again concavity assumptions. Namely that concavity on r is strong enough to 
compensate for concavity in C 

- Eli\x\ 2 < r(x) - sup C(x,y)p M < r(x) < k' -K'^x] 2 , (9.4) 
y 

- 2^; < D 2 r(x) - sup (D 2 C(x,y)) +PM < D 2 r(x) + sup (D 2 C(x,y)) _p M < -2K[, (9.5) 

y 1 y 

as symmetric matrices, where the positive and negative parts are taken componentwise. As for 
regularity, we will use 

D 3 r- sup (D 3 C(-,y)) p M , D 3 r + sup {D 3 C(; y))_p M € L°°(M d ). (9.6) 
y y 

The initial data is still supposed to concentrate at a point x° following (|3.8p - ([3.1ip . But 
because persistence, i.e. that n £ does not vanish, is more complicated to control, we need two 
new conditions 

r(x°) > 0, (9.7) 



n e (t,x)dx < p° M - (9.8) 
We also need a compatibility condition with R 

41? < K[ < < ALj . (9.9) 

The interpretation of our assumptions is that the intra-specific growth rate r dominates 
strongly the competition kernel. This avoids the branching patterns that are usual in this kind 
of models [22], [TTJ, [251 ES] • O ur concavity assumptions also implies that there is no continuous so- 
lution iV to the steady state equation without mutations N(x) (r(x) — J Rd C(x, y)N(y) dy) = 0. 
This makes a difference with the entropy method used in }18] as well as the positivity condition 
on the kernel that, compared to (|9.3p . also involves r(x). 

Our goal is to prove the following results 
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Theorem 9.1 (Convergence) With the assumptions ([Q ]) - ([93|l and (pT8]h (|978|) — ([979]) . the 

solution n e to (19. If) satisfies, 

0<Pe(t):= / n e (t,x)dx < p M a.e. (9.10) 

and £/iere is a subsequence such that 

pJt) — ^ pit) in weak-* L°°(R + ), < pit) < p M a.e. (9.11) 

Furthermore, we have weakly in the sense of measures for a subsequence n t 

anc? i/ie pair (x(i),/o(t)) a/so satisfies 

p(t) [r(x(t)) - p(t)C(x(t), x(t))] > 0. (9.13) 

With the assumptions of Theorem 19.11 we do not know if p e converges strongly because we do 
not have the equivalent of the BV quantity in Theorem 13. II We can only prove it with stronger 
assumptions. This is stated in the 

Theorem 9.2 (Form of canonical equation) We assume (|3.8p ~ (|3.1ip and f)9.3[) — fj9.9[) . Then, 
for the function C 2 -function u(t,x) given below in (I9.24D with D^.u G L^ C (M + ; L OD (R d )), x 6 
W l '°°{M. + ) satisfies 

x(t) = (-D^xit)))' 1 • [V x r(x(t)) -p(t)V x C(x(t),x(t))] , (9.14) 
with initial data x° given in (|3.1ip . Furthermore, p e converges strongly and we have p(t) E 

iy 1 ' oo (R+), 

r(x(t)) - p(t)C(x(t),x(t)) = 0, (9.15) 
r(x(t)) > r{x°)e~ Kt , p{t) > p°e' Kt . (9.16) 

We may find some kind of gradient flow structure for the canonical equation when C(x,y) is 
symmetric and obtain 

Theorem 9.3 (Long time behavior) We make the assumptions of Theorem C(x,y) = 
C(y,x) and 

x I—)- &(x) := lnr(x) — \nC(x,x) is strictly concave in the set {r > 0}. (9-17) 
Then, as t — )• oo, p(t) — > p^ > 0, x(t) — > Xqo and Xoo is the maximum point of 
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9.1 A-priori bounds on p e 



9.1 A-priori bounds on p e 

The main new difficulty with the competition model comes from a priori bounds on the total 
population. In particular it is not known if there are BV quantities proving that p e {t) converges 
strongly. Even non extinction is not longer automatic. 

One side of the inequality (|9.10p is given by n e > 0, for the other side we integrate (|9.ip over 
R d and use JOD to write 



— / n,(t,x)dx = - [ n<(t,x)r(x)<fc-i / / n,(t,x)C(x,y)n,{t, y) dydx 
dt ./rod e ./rod e 



> - / n e (t, x)r(x) dx 



J^dn e (t,x) dx 
PM 



e 

therefore (and even though r can change sign) we conclude thanks to 

n e (t, x) dx < pm- 

9.2 Passing to the limit 

The proofs of the remaining parts of Theorems are close to those already written before. We 
only give the main differences here. They rely again on the WKB ansatz u t := eln(n e ). We 
obtain as before that u t satisfies the Hamilton-Jacobi equation 



(9.18) 



j §- t u e (t,x) = r(x) - / C(x,y)n € (t,y)dy + \Vu e \ 2 + eAu € , t > 0, x G R d , 
[ Ue (t = 0) = u° e . 

Similarly to Lemma 15. II we can prove the 

Lemma 9.4 With the assumptions of Theorem \9.1\. we have for all t > 

-La ~ Li\A 2 ~ e2dL.it < u e (t,x) < L -L~i\x\ 2 + (l{' + 2de~L^j t, 

-2L_i < D 2 u e (t,x) < -2Li. 

Proof. The first line holds because the right (resp. left) hand side of the inequality is a super 
(resp. sub) solution thanks to assumption (|9.4p and using the control of n £ by pu- For the 
second line, the upper and lower bound use the maximum principle on the equation for D 2 u £ 
and the compatibility conditions (|9.9p as in section 15. 21 □ 

From the regularity obtained, it follows that we can extract a subsequence such that u e (t, x) — > 

E— >0 

u(t,x), locally uniformly as in section [5\3l We also obtain from Lemma 19.41 

-Lo -IaM 2 < u(t,x) < L - Ii|x| 2 + K' t, -2L X < D 2 u(t,x) < -2L~i a.e. (9.19) 

and that u satisfies, in the viscosity sense (modified as in [4l 124^13]). the equation 

f §i u = r ( x ) ~ P(t)C{x, x(t)) + | Vn| 2 , ^ 2Q 

u(t,x) < 0. 
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The constraint is now relaxed to an inequality because we know that the total mass is bounded 
but we do not control the mass from below at this stage. In the framework of Theorem 19. 2\ we 
prove later on that the constraint is always an equality (see (|9.24p ). 
It might be that p(t) vanishes and then x(t) does not matter here, nevertheless we still have 



max 



u(t,x) = u(t,x(t)). (9.21) 



Using the control (|9.19p and this Hamilton-Jacobi equation we obtain (19.120 with the same 
arguments as in Section 15.31 The new difficulty is that p(t) might vanish in particular when the 
constraint is strict max R d u(t, x) < 0, an option that we will discard later. Because of that, we 
also obtain the restriction on times in (|9.13p which can be completed as (in the viscosity sense) 

j t u(t,x(t)) =r[x(t)) - p(t)C(x(t),x(t)). (9.22) 



We also have, as in section 19.11 

4*W=/ MMH^-f / n € {t,x)C(x,y)n € {t, y) dydx. 
Passing to the weak limit (integration by parts and using boundedness of p t ), we find that 

p(t)r (x(t)) = w-lim / / n e (t,x)C(x,y)n e (t,y) dydx > p{t) 2 C{x{t),x{t)) . (9.23) 
This proves (|9.13f) and concludes the proof of Theorem 19.11 



9.3 Form of the canonical equation 

We continue with the proof of Theorem 19.21 and we begin with the derivation of (|9.14p . 

The third derivative of u e is bounded using assumption (|9.6p and following the same arguments 
in Section 16.11 Then similarly to Section 16.21 we have the regularity D%u 6 L°°((0,T) x R d ) for 
all T > 0, g u and Df xx u £ L^ C (R + x R d ). 

The canonical equation (19. 141) is established a.e. as in section 16.21 using the maximum points 
of u € and passing to the limit. From (|9.10p . ()9.14j) and (|9.19|) . we next obtain that \-^x(t)\ is 
uniformly bounded. 



9.4 Persistence 

Now we prove that u(t,x(t)) = for all t > and that p(t) > a.e. t. We cannot obtain 
this directly and thus we begin with proving r{x{t)) > 0. We prove indeed the first part of the 
inequality (|9. 16j) . 



We prove this by contradiction. We suppose that to is the first point such that r(x(to)) = 0. We 
notice that x(t) being lipschitzian and using assumption (|9.7p . for all t < to, we have r(x(t)) > 0. 
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9.4 Persistence 



Therefore with assumption (|9.4p we deduce that x(t) remains bounded for t £ [0, to]. Using (|9. 14[) 
and (I9.19P we have 

at 

= V x r(x(t)) ■ {-D 2 u{t,x{t)))' 1 ■ [V x r(x(t)) - p(t)V x C(x(t),x(t))] 
>-p(t)\V x r(x(t))\\V x C(x(t),x(t))\. 
Consequently, using (|9.13p . we obtain 

ir{x(t)) > -r(m)\^r(m\- VxC{m,m)l 



dt y y " ~ K v x v v C(x(t),x(t)) 

Moreover we know that x(t) remains bounded for t £ [0, to] an d thus we have inf C(x(t),x(t)) > r]2 > 0. 

te[o,to] 

We conclude that, for K a positive constant, 

■^r(x(t)) > -Kr(x(t)), for < t < t . 

Starting with r(x°) > according to (|9.7|) . this inequality is in contradiction with r(x(to)) = 0. 
Therefore for all t > we have r(x(t)) > and thus this inequality is true for allt > 0. Thereby 
we obtain the first part of f|9. 16|) . From the latter and using again (|9.4p . we also deduce that 
x(t) remains bounded for all t > 0. 

Next we use (|9.13p . (|9.22p and the positivity of r(x(t)) to obtain 

u(t,x(t))-v°(x°) = J r(x(t))-p(t)C(x(t),x(t))dt > l- p[t)=0 r(x(t)) > 0. 

We deduce, using (|3. 11 [) . that u(t,x(t)) = for all t > 0. Thus the equation on u is in fact 

f &u = r{x)-p(t)C(x,x(t)) + |Vn| 2 , 

< (9.24) 
[ max Rt j u(t, x) = 0. 

This identity also proves the identity (|9.15p and thus that (|9.23p holds as an equality, which 
is equivalent to say that the weak limit of p e {t) is in fact a strong limit. 

We may now use (19151) to conclude that pit) is also bounded in To do so, we first 

differentiate fl9. 151) and find again some kind of gradient flow structure 

Vr (x(t)) • x(t) - p(t)C(x(t), x(t)) - pit) [V x C[x(t), x(t)) + V v C{x(t), x(t))] x(t) = 0. 

With (HHD, it follows that 

x(t) ■ (~D 2 u) ■ x(t) = p(t)C(x(t),x(t)) + p(t)V y C(x(t),x{t))x(t). (9.25) 

From the uniform bounds proved before, there is a constant 773 such that 

V v C(xit),x(t)) -xit) 
V Ll ' \ L < m- (9-26) 
C{xit),x{t)) 

Using the latter and (|9.19p we conclude that, for K a positive constant, 

Pit) > -Kpit). 

Thus we obtain (|9.16|) . This completes the proof of Theorem 
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REFERENCES 



9.5 Long time behavior 

It remains to prove Theorem 19.31 Assuming that C(x,y) is symmetric, we can find a quantity 
which is non-decreasing in time, which replaces the quantity I in section [6T3l We compute, from 
the relation ([9T25|> . 

j t [p\t)c{x(t),x(t))] mt))- 1 = p(t)c[x(t),m)+p(t)(y y c)m 

= x(t) ■ (~D 2 u) ■ x(t) > 0. (9.27) 

As t tends to infinity, we may consider a subsequence t n such that p(t n ) — > p^, x(t n ) — > Xoo. 
From ()9.27p . we may also assume x(t n ) — > 0. Therefore 

Vr(x oo ) PoO^xC^Xqo, Xqo), ^(^Oo) PooC^Xqq, Xqq). 

From these relations, we first conclude that poo is positive because r is concave and its gradient 
only vanishes at a point where r is positive. 
Then we combine the relations and conclude that 

Vr(xoo) V x C(x 

OO) 3-oo ) 
'"(■Eoo) CX^ooj^oo) 

The assumption (I9.17P then concludes on the uniqueness of such a point x^ and thus on the 
convergence of x(t). □ 
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